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ABSTRACT 

Relatively  little  attention  has  been  paid  to  the  potential  of  polarization  techniques  to  provide  additional  information 
for  the  mapping  of  earth  resources,  compared  to  the  published  work  concerning  other  active  and  passive  i emote 
sensing  systems.  Recently,  a  substantia!  number  of  polarized  light  images  of  a  variety  of  terrestrial  scenes  have 
been  obtained  from  the  Space  Shuttle.  A  boresighted  pair  of  Hasselblad  cameras  was  used,  in  which  polarization 
filters  were  fitted.  The  polarization  directions  were  perpendicular  to  each  other  for  the  two  cameras.  Image  pairs 
were  acquired  with  one  image  being  of  maximum  intensity,  and  the  other  showing  minimum  intensity.  Selected 
pairs  of  images,  obtained  using  black  and  white  films,  were  digitized.  The  images  were  registered  and  compared, 
using  digital  image  analysis  techniques.  Differences  due  to  polarization  were  observed,  these  included  intensity  and 
contrast  differences,  together  with  differences  in  the  spatial  frequency,  orientation  and  population  of  observable 
contrast  boundaries.  It  was  found  empirically  that  some  digital  image  analysis  techniques  enhanced  the  differences. 
A  comparison  of  the  enhanced  difference  images,  obtained  from  the  polarization  pairs,  with  the  intensity  data 
that  would  have  been  obtained  without  the  aid  of  polarization  filters,  confirms  that  a  considerable  degree  of  new, 
useful  infoimation  may  be  obtained  by  the  use  of  polarization  techniques.  This  new  information  was  most  helpful 
in  better  defining  the  observed  ground  features  in  the  cases  studied.  This  is  because  the  better  definition  of  contrast 
boundaries,  improvement  of  contrast  across  boundaries,  improvement  of  shadow  detail  and  reduction  of  noise 
level  increases  the  useful  information  in  an  image,  and  improves  its  interpretability. 

1.  INTRODUCTION 

Although  many  laboratory  studies  of  polarization  have  been  made  compared  to  research  studies  in  the  unpolarized 
visible  and  infrared  regions  there  has  been  little  activity  to  attempt  to  translate  polarization  phenomena  into  a  technique 
for  remotely  mapping  earth  resources  or  to  study  meteorological  or  oceanographic  phenomena.  Egan1  deduced 
from  laboratory  measurements  on  vegetation  that  information  gained  from  polarization  measurements  could  assist 
in  the  mapping  of  agricultural  crops.  Vanderbilt,  Grant  and  Daughtry2,  Egan’,  Slater4,  and  Talmadge  and  Curran5 
point  out  the  potential  value  of  polarization  techniques  in  remote  sensing.  Coulson6  and  Coulson,  Bouricius  and 
Gray"  have  made  ground  and  laboratory  measurements  of  the  view  angle  and  illumination  angle  dependence  of 
t’  polarization  of  light  reflected  from  surfaces.  More  recent  measurements  have  been  made  by  Vanderbilt,  Grant 
and  Daughtry2.  Halajian  and  HallockK  have  investigated  polarimetric  mapping  as  a  remote  sensing  tool.  Curran9 
has  reported  on  the  use  of  aerial  polarimetric  photographic  methods  in  soil  moisture  investigations. 

More  recently,  the  advent  of  space  shuttle  has  enabled  the  exploratory  use  of  polarized  light  photography  of  the 
earth  for  the  preliminary  assessment  of  the  utility  of  polarization  methods  for  earth  resources  mapping,  as  reported 
by  Coulson,  Whitehead  and  Campbell1".  A  recent  book  by  Coulson11  provides  a  thorough  coverage  of  present 
knowledge  regarding  polarization  effects  of  the  atmosphere. 

In  this  paper,  we  report  upon  fm-ther  preliminary  investigations  of  polarized  light  photography  of  the  earth  for 
earth  resources,  meteorological  and  oceanographic  observation.  Due  to  the  constraints  of  time  and  space,  only 
part  of  our  investigations  are  reported  here. 
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2.  DISC  IJSSION  OF  PAST  WORK 

The  mean  polarization  state  of  radianee  reflected  from  a  given  feature  may  change  with  time,  due  to  the  change 
with  time  in  the  orientation  of  the  mean  electric  vector.  The  polarized  wave  is  the  resultant  of  the  interaction  of 
the  ensemble  of  waves  propagated  in  the  same  direction  but  w  ith  different  phase  and  amplitude  at  different  orientations 
about  the  direction  of  propagation.  In  the  case  of  linear  polarization,  the  orientation  of  the  plane  of  the  resulting 
polarized  wave  about  the  direction  of  propagation  does  not  vary  with  time.  In  the  case  of  circular  polarization, 
the  ensemble  of  propagated  waves  produces  the  same  effect  as  two  orthogonal  waves  of  equal  amplitude  which 
are  out  of  phase  by  a  constant  phase  angle.  The  resulting  wave  is  of  constant  amplitude  but  the  plane  of  the  electric 
vector  rotates  with  time  about  the  direction  of  propagation.  In  the  case  of  elliptical  polarization,  the  amplitude 
of  the  resultant  wave  v  aries  with  time  as  the  plane  of  the  resultant  wave  rotates  about  the  direction  of  propagation. 
These  three  states  are  shown  in  Fig.  1.  Here,  we  are  looking  down  the  axis  of  propagation  at  the  locus  traced  out 
by  the  amplitude  of  the  resultant  electric  wave.  The  orientation  of  the  semimajor  ellipse  of  the  polarization  vector 
w  ill  depend  upon  the  magnitude  and  polarity  of  the  two  components  of  the  resultant  electric  field  vectors  into  which 
the  ensemble  of  wave,  propagating  in  the  same  direction  can  be  decomposed. 

In  order  to  define  the  state  of  polarization,  we  must  consider  those  properties  that  we  need  to  measure  that  will 
characterize  its  polarization.  G.  G.  Stokes,  in  1852,  introduced  four  quantities  that  are  functions  only  of  the  observables 
of  the  electromagnetic  wave  and  are  known  as  Stokes  parameters.  A  useful  description  is  in  Hecht and  in  Egan  T 

Egan*  gives  a  definition  of  the  four  Stokes  parameters  in  terms  of  the  amplitudes  of  the  electromagnetic  waves 
in  the  x  and  \  directions  which  are  mutually  orthogonal,  with  .v  being  in  the  plane  of  the  electric  vector  for  minimum 
polarization  and  y  being  in  the  plane  of  the  electric  vector  for  maximum  polarization,  which  corresponds  to  the 
directions  of  maximum  and  minimum  intensity,  respectively,  as  an  analyzer  is  rotated  in  the  beam. 


/  = 

<A2  +  A2>  =  A2 

(1) 

Q  - 

<A2-A2> 

(2) 

u  = 

<2  A  xA  ,.  cos  y  > 

(3) 

v  = 

<2  A  xA  v  sin  y  > 

(4) 

This  is  a  time  average  of  the  definition  given  by  Chandrasekhar Here,  <  >  indicates  time  averaging.  A  v  and 

•1 1  are  the  magnitudes  of  the  representative  average  electromagnetic  waves  in  the  .v  and  y  directions,  A  2  is  intensity 
and  y  is  the  phase  angle  between  A  (  and  A  v.  The  radiant  intensity  /  is  related  to  the  othr .  three  Stokes  parameters 
by  the  equation 

r-  =  q:  +  u2  +  v2 .  (5) 

Q  is  the  difference  in  radiant  intensity  between  the  mutually  orthogonal  .v  and  y  directions  used  to  specify  A  x  and 
•I,.  In  the  case  where  the  .v  and  y  directions  represent  the  directions  of  the  electric  vector  for  maximum  and 
minimum  intensity,  respectively,  then 


^ma\  ^min 


whereas,  for  this  case 


/ 


I  +f 

max  '  min 


(6) 

(7) 


In  the  case  of  partially  polarized,  linear,  quasimonochromatic  light,  if  we  rotate  an  analyzer  in  the  beam,  then 
there  will  be  an  orientation  at  which  the  transmitted  irradiancc  is  a  maximum.  This  will  be  at  right  angles  to  the 
orientation  where  it  is  at  a  minimum.  In  this  cac  \  the  maximum  intensity  will  consist  of  the  sum  of  the  polarized 
(lp)  and  unpolarized  (/„)  components,  while  the  minimum  component  will  consist  only  of  the  unpolarized 
component  I u. 
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Thus  0  =  (/„  +  lu)~  (/„) 


(8) 


whereas 


=  / 


max 


-  / 


min  ’ 


(9) 


1  ~  Aiax  +  ^min  ~  lp  +  2A< 


Then  the  degree  of  polarization  (P)  of  the  beam  is  given  by 


</4:  -  Al> 


P(* o)  =-■  100  ,v  =  100  t  =  100  * 


2 

/ 


,  + 

The  shape  of  the  polarization  ellipse  is  given  by 

.  ,  b 

x  =  tan  1  -  , 
a 

where  a  and  b  are  the  semimajor  and  the  semiminor  axes  respectively.  Also 

_  _1_  tan  1  V 
X  ~  2  (Q2  +  U2)'  2  ’ 


Anax _ An  in 

^max  Alin 


(10) 


(11) 


(12) 


(13) 


W  hether  the  polarization  is  right-handed  or  left-handed  depends  upon  the  sign  of  V.  The  orientation  of  the  polarization 
ellipse  with  respect  to  a  reference  direction  is  given  by 


(14) 


Methods  have  been  summarized  by  Talmadge  and  Curran5  for  measuring  reflected,  polarized  light  with  linear 
polarizing  filters  set  at  different  orientations  to  determine  the  Stokes  parameters.  As  observed  by  Coulson,  Whitehead 
and  Campbell1"  and  by  Talmadge  and  Curran5,  the  Stokes  parameter  Kis  generally  very  small  after  reflection  from 
the  earth's  surface.  V  represents  the  amount  of  circularly  polarized  light  (e.g.,  Egan-1).  Indeed,  Spohr14  reports 
that  the  contribution  of  the  circular  polarized  component  to  the  total  intensity  is  tw^  orders  of  magnitude  less  than 
the  contribution  of  the  linear  polarized  component.  Since  V  is  very  small,  Prosh,  Hennings  and  Raschke15 
determine  the  Stokes  parameters  /,  Q ,  and  U  from  radiance  measurements  using  polarizers  in  front  of  separate 
radiometers.  They  compute  the  components  of  the  Stokes  vector,  where  the  polarization  filters  are  at  respective 
orientations  about  the  optic  axis  of  the  detectors  of  0°,  60°  and  120°.  The  Stokes  vector  components,  expressed 
in  terms  of  the  transmitted  radiance  In,  /Wl,  are 


1  =  \  (Ai  +  +  hm) 

Q  ~  J  (2 A)  _  4)  “  A:,i) 

k  ~  ( A  20  “  Ao) 


(15) 

(16) 

(17) 


Prosh,  Hennings  and  Rashke  converted  the  output  voltages  from  the  three-detector  video  polarimeter  to  color 
coordinates  for  display  using  a  color  video  monitor.  The  azimuth  of  polarization  (x)  defines  tint,  and  the  degree 
of  polarization  (P)  defines  hue.  They  reported  that  their  instrument  and  technique  provided  useful  quick-look  data 
and.  when  used  on  an  airborne  platform,  enhanced  the  discriminability  of  features  on  the  ground.  Walraven16  used 
a  35mm  camera  to  obtain  sets  of  four  slides  of  scenes  using  a  polarizing  filter  at  0°,  45°,  90°,  and  135°  relative 
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orientation  of  the  polarization  direetion  about  the  optic  axis  of  the  camera.  The  intensity  of  the  radiation  that 
passes  through  the  polarization  filter  is  given  by 


/'  (0)  =  1/2 (/  +  Q  cos  20  +  U  sin  20)  ,  (18) 

where  0  is  the  angle  of  the  transmission  axis  of  the  polarizer  about  the  optic  axis  of  the  camera,  with  respect  to 
the  horizontal.  If  measurements  are  successively  made  with  0  =  0°,  45°,  90°,  and  135°,  then  the  resulting  intensity  is 


/'( 0)  =  l/2(/  +  Q)  , 

(19) 

/' (45)  =  !/2(/  +  U)  , 

(20) 

7’ (90)  =  1/2(7  -  O)  , 

(21) 

/'  ( 1 35)  =  1/2(7  -  U)  . 

(22) 

From  these  equations,  Walraven  deduced  the  Stokes  parameters: 

7  =  1/2  [/'((»  +  7' (45)  +  7' (90)  +  7'(135)] 

(23) 

Q  =  7'( 0)  -  7'  (90) 

(24) 

U  =  7' (45)  -  7  (135) 

(25) 

With  the  polarizer  set  at  each  of  the  four  angles  defined  above,  the  four  photographic  images  were  obtained  with 
the  same  f-stop.  They  were  then  digitized  and  registered.  Images  of  7,  P,  £  were  obtained: 

P  =  j  (Q1  +  U2)'2  (26) 

The  Stokes  parameter  V  is  very  small,  and  the  azimuth  of  polarization  with  respect  to  the  horizontal  (£)  is  given 
by  equation  (14). 

Walraven  showed  images  of  7,  (1-P)  and  £  for  four  natural  scenes.  He  used  (1-P)  instead  of  the  degree  of  polarization 
P,  since  the  degree  of  polarization  appears  to  be  inversely  proportional  to  the  albedo.  He  showed  that  the  polariza¬ 
tion  images  provided  new  and  useful  information  about  a  scene  beyond  that  available  from  an  intensity  image. 

As  Walraven 16  points  out,  the  measurements  made  in  the  laboratory  are  considerably  more  complicated  when  made 
in  the  field.  Illumination  is  no  longer  from  a  point  source,  but  partly  from  scattered,  diffuse  sky  illumination  which 
is  polarized  to  an  extent  which  depends  on  the  angular  location  of  each  elemental  source  of  irradiance  in  the  sky. 
The  irradiance  from  the  sun  is  unpolarized,  although  that  from  the  sky  is  polarized  to  an  extent  which  depends 
on  sun  angle  and  upon  the  reflectance  of  terrain,  as  shown  by  Coulson."'7  Variations  in  composition, 
heterogeneity  and  topography  may  also  arise  at  both  a  resolvable  and  at  a  subresolution  level  (e.g.,  Duggin18). 
Walraven  also  points  out  that  since  man-made  surfaces  are  unusually  smooth,  they  exhibit  higher  than  usual 
polarization.  Darker  surfaces  tend  to  exhibit  larger  polarizations  because  the  single  scattering  component  is  larger 
than  that  for  light  surfaces,  where  multiple  scattering  prevails.16 

Coulson,  Whitehead  and  Campbell10  have  observed  that  there  is  little  information  regarding  orbital  observations 
of  polarization  of  light  reflected  from  the  earth’s  surface,  although  there  is  some  data  on  the  reflection  of  polarized 
light  from  the  earth’s  atmosphere. Iv  Coulson,  Whitehead  and  Campbell10  studied  polarization  measured  from 
digitized,  registered  polarization  photographic  pairs  obtained  from  space  shuttle  flight  STS-5 1 .  These  authors  analyzed 
image  data  obtained  from  a  synchronized  hand-held  camera  pair  in  which  both  cameras  were  attached  to  a  common 
frame  and  bore-sighted.  The  linear  polarization  filters  had  polarization  axes  oriented  at  right  angles  to  each  other 
about  the  optic  axis.  The  cameras  were  pointed  at  terrestrial  targets  and  rotated  about  the  optic  axis  until  one  camera 
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showed  maximum  intensity,  and  the  other  showed  minimum  intensity  through  the  viewfinders:  such  judgments 
were  empirical.  The  photographs  were  subsequently  digitized  and  analyzed  to  provide  polarization  images.  The 
authors  concluded  that  there  was  potential  for  the  application  of  polarization  photography  from  orbit  in  studies 
of  the  earth.  They  observed  degrees  of  polarization  of  up  to  0.5.  The  degree  of  polarization  varied  systematically 
across  the  image.  It  was  reported  that  the  polarization  image  made  some  patterns  discernable  which  were  not  discernible 
on  the  intensity  images.  Further,  it  was  reported  that  the  use  of  polarization  improved  the  signal-tc- 
noise  ratio. 

In  this  paper,  we  report  an  analysis  of  the  black-and-white  film  image  pairs  from  space  shuttle  mission  STS-51A. 

Since  image  pairs,  representing  maximum  and  minimum  intensity  only,  were  obtained,  it  is  not  possible  to  calculate 
all  of  the  Stokes  parameters.  If  we  assume,  on  the  basis  of  prior  work101415  that  the  circular  component  of 
polarization  of  radiance  reflected  from  the  earth  is  close  to  zero,  then  we  could  obtain  more  complete  information  if 
we  have  sequences  of  three  images  with  the  polarization  filter  oriented  at  three  angular  positions  about  the  optical  axis 
of  the  camera  system.  In  addition  to  the  polarization  filters  whose  polarization  directions  are  orthogonal,  an  additional 
channel  with  the  polarization  filter  oriented  at  (say)  45°  to  the  maximum  could  be  used. 

Howes  er,  the  data  analyzed  were  collected  for  exploratory  purposes.  The  analyses  of  this  two-channel  data  clearly 
show  that  as  found  by  Coulson,  Whitehead  and  Campbell10  there  is  ample  reason  to  expect  a  great  deal  of  benefit 
for  earth  resources  investigations  from  the  orbital  acquisition  and  analysis  of  polarized  light  imagery  for  analysis 
m  conjunction  with  earth  resources  studies. 

3.  EXPERIMENTAL  DATA 

fhe  polarization  image  pairs  were  obtained  on  space  shuttle  mission  STS-51A  (October  1984).  Photography  was 
obtained  using  a  boresighted  pair  of  Hasselblad  cameras,  which  were  synchronized  to  produce  simultaneous  exposures. 
Polarization  filler s  were  fitted  to  both  cameras,  and  the  directions  of  polarization  were  orthogonal  to  each  other 
about  the  optic  axis  of  the  camera  system.  The  camera  pair  was  rotated  about  its  optic  axis  by  the  astronaut 
photographer  so  that  when  a  target  of  interest  was  selected,  the  image  through  one  viewfinder  was  of  maximum 
brightness,  and  the  image  through  the  other  viewfinder  was  of  minimum  brightness.  Thus,  two  images  were  obtained: 
one  of  maximum  intensity  and  the  other  of  minimum  intensity.  We  used  second-generation  film  that  was  processed 
to  positive  transparency.  There  is  no  way  of  knowing  the  precision  in  selecting  the  maximum  and  minimum  intensity, 
nor  of  knowing  the  change  in  polarization  caused  by  the  space  shuttle  windows.  Therefore,  the  selection  of  angles 
w  as  assumed  to  be  exactly  correct,  and  the  potential  impact  of  the  space  shuttle  windows  on  polarization  was  ignored, 
fhe  images  were  digitized  using  a  Diconix  digitizer  and  VAX  11/780  computer  system  at  NASA’s  Johnson  Space 
(  enter.  Unfortunately,  the  film  had  no  sensitometry  strip;  thus  the  films  from  both  cameras  were  assumed  to  be 
identically  processed,  so  that  the  calibration  of  focal  plane  irradiance  to  photographic  density  would  be  the  same  in  each 
case,  f  ollowing  this  assumption,  the  degree  of  polarization  could  be  obtained  from  digitized,  registered,  image  pairs. 
W  e  used  8-bit  quantization  during  image  digitization  and  registered  the  images.  The  digitized,  registered,  image 
pairs  were  than  manipulated  using  standard  image  analysis  techniques  as  described  in  the  following  section. 

4.  ANALYSIS  AND  RESULTS 

We  selected  20  images  from  the  black  and  white  image  pairs  obtained  on  Space  Shuttle  Mission  STS-51 A  for  further 
-tudy.  The  digitized  images  were  registered  using  nearest  neighbor  interpolation,  so  the  radiometry  of  the  images 
would  not  be  modified.  While  the  different  images  obtained  different  radiometric  information,  it  was  considered 
that  spatial  information  in  the  form  of  contrast  boundaries  and  texture  would  be  different  in  the  two  images.  While 
only  the  degree  of  polat  ization  (i.e.,  the  first  two  Stokes  parameters)  may  be  obtained  from  an  image  pair,  information 
on  contract  boundaries  and  on  texture  differences  between  components  of  the  image  pair  should  be  able  to  be  extracted, 
given  the  appropriate  analysis.  To  this  end,  after  carefully  digitizing  the  image  pairs  using  the  1CON1X  facility 
at  NASA  Johnson  Space  Center,  we  performed  digital  image  analysis  on  the  1DIMS  installation  at  NASA  Ames 
Research  Center;  on  the  F.RDAS  installation  at  State  University  of  New  York  ESF,  Syracuse;  and  on  the  FS  System 
600  at  the  Naval  Ocean  Research  and  Development  Activity,  Stennis  Space  Center,  Mississippi. 
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Each  image  pair  was  scaled  and  individually  displayed  together  with  images  of  the  sum  and  the  difference;  the 
ratio  and  difference  divided  by  the  sum.  In  an  empirical  investigation,  we  also  combined  the  individual  images, 
the  sum,  difference,  and  ratio,  and  difference  divided  by  the  sum  images  to  provide  input  to  a  principal  component 
analysis.  Of  course,  the  input  images  were  strongly  correlated,  but  the  principal  component  images  did,  in  some 
cases,  provide  a  useful  data  enhancement.  We  also  obtained  beneficial  pictorial  image  information  by  performing 
a  principal  components  analysis,  performing  histogram  equalization  on  the  principal  components  images,  and  then 
applying  the  inverse  transformation  matrix.  The  results  we  obtained  are  extremely  valuable,  but  we  cannot  completely 
explain  at  this  time  why  such  analyses  are  so  effective. 

Further  work  consisted  of  calculating  the  Fourier  spectra  of  each  of  the  images  and  of  calculating  the  difference 
in  Fourier  spectra  between  the  components  of  the  image  pair.  The  difference  in  Fourier  spectra  will  show  the 
difference  in  both  contrast  boundaries  and  edge  features  provided  by  contrast  differences  between  the  two  images. 
A  significant  azimuthal  anisotropy  was  observed  to  occur  in  the  spatial  frequency  plane.  The  appearance  of  this 
property  was  not  altogether  surprising,  since  differences  in  the  shadow  components  between  the  polarization  image 
pair  will  occur  due  to  the  polarization  characteristics  of  the  ground  and  the  polarization  characteristics  of  the 
atmosphere,  and  will  be  strongly  affected  by  sun-target-sensor  geometry.  Shadowed  areas  will  be  almost  entirely 
diffusely  illuminated,  so  that  differences  in  the  depth  of  shadow'  will  be  recorded  as  difference  in  the  degree  of 
polarization.  Indeed,  one  objective  of  our  investigation  is  to  determine  whether  and  to  what  extent  polarized  light 
may  be  used  to  enhance  detail  in  shadowed  areas.  Thus,  the  spatial  frequency  information  contained  in  each  image 
of  the  polarization  pair  will  be  impacted  by  the  topography  and  optical  characteristics  of  the  ground,  and  by 
atmospheric  conditions,  as  well  as  by  view  and  illumination  geometry.  Further,  water  bodies  will  strongly  polarize 
reflected  light;  the  degree  of  polarization  will  depend  upon  the  azimuth  angle  between  the  planes  containing  the 
detected  ray  and  the  corresponding  radiance  from  the  sun.  Appropriate  selection  of  illumination  and  imaging 
conditions  may  therefore  reasonably  be  expected  to  yield  useful  hydrological  information.  Oblique  views  of  anisotropic 
ground  scenes,  illuminated  from  the  side  will  generally  result  in  a  greater  degree  of  difference  between  the  Fourier 
spectra  of  polarization  image  pairs. 

An  example  of  an  image  obtained  over  a  land  area  containing  hydrologic  features  (Niger  River  in  Mali)  is  shown 
in  Fig.  2.  The  components  of  the  polarization  pair,  images  A  and  B,  the  difference  image  (A-B),  the  sum  image  (A  +  B), 
the  ratio  image  ( A/B)  and  the  degree  of  polarization  image  (A-B)/(A  +  B)  are  shown.  The  Niger  River  is  displayed  well 
in  Fig.  2a,  but  is  almost  invisible  in  Fig.  2b,  the  other  component  of  the  polarization  pair.  The  difference  image 
(Fig.  2c)  provides  slight  enhancement  of  the  river  system,  but  while  the  sum  image  (Fig.  2d)  is  better  than  image  B, 
it  has  less  hydrological  information  than  does  image  A,  as  one  might  expect.  The  degree  of  polarization  image  (Fig.  2e) 
does  not  enhance  the  river  system,  but  does  enhance  sand  dune  structure.  The  ratio  image  (Fig.  20  does  not  provide 
more  information  on  the  river  than  the  degree  of  polarization  image. 

A  principal  component  analysis  was  performed  with  input  consisting  of  the  individual  images,  the  sum,  difference,  the 
ratio,  and  polarization  images.  The  first  and  second  principal  components  images  are  shown  in  Figs.  2(g)  and  2(h). 
Clearly,  much  more  complete  information  on  the  Niger  River  and  its  tributaries  is  contained  in  the  second  principal 
component  image  than  in  any  of  the  input  images  (Figs.  2(a)-2(f)).  A  further  procedure  was  also  investigated,  since 
the  input  images  were  highly  correlated.  The  principal  components  were  histogram-equalized,  and  the  inverse 
transformation  matrix  was  applied.  The  resulting  first  inverse  transform  image  (essentially  an  enhanced  version 
of  image  A,  containing  information  from  all  of  the  other  input  bands)  is  shown  in  Fig.  2(i).  It  is  not  clear  why 
additional  information  was  obtained  from  this  treatment  of  the  strongly  correlated  input  images.  Flowever,  we  present 
this  empirical  data  as  evidence  that  digital  image  analysis  can  provide  a  means  of  extracting  maximum  information 
from  the  differences  between  the  input  polarization  image  pair.  We  then  present  this  information  in  a  single  image 
that  will  facilitate  mapping  a  specific  feature  of  interest. 

We  also  explored  the  differences  in  the  spatial  frequency  spectra  of  features  delineated  by  contract  boundaries  in 
the  two  images.  The  maximum  and  minimum  intensity  images,  which  correspond  to  the  minimum  and  maximum 
polarization,  respectively,  contain  contrast  and  contrast  boundary  information  that  is  different  in  each  image.  The 
delineation  and  recognition  of  features  will  depend  largely  on  the  recognition  of  contrast  boundaries.  If  the 
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intensity  information  differs  from  one  image  to  another,  then  the  contrast  boundary  information  will  differ  between 
the  images.  Once  the  contrast  across  a  boundary  falls  below  a  certain  level,  the  boundary  may  not  be  detected. 
The  subtraction  of  the  Fourier  spectrum  of  one  image  from  that  of  the  other  image  constituting  the  polarization 
pair  will  show  the  difference  in  the  spatial  frequencies  at  which  contrast  boundaries  occur  in  the  image  pair,  and 
will  show  the  azimuthal  distribution  with  which  they  occur.  The  frequency  of  occurrence  of  boundaries  appearing 
with  given  spatial  frequencies  at  given  azimuthal  orientations  will  be  shown  by  the  intensity  of  points  represented 
in  the  frequency  plane.  The  difference  in  the  Fourier  transforms  of  the  two  images  is  shown  in  Fig.  2(j).  Clearly, 
there  is  not  only  a  difference  in  the  frequency  of  occurrence  in  contrast  boundaries  between  the  images,  but  there 
is  also  an  azimuthal  anisotropy.  One  might  expect  this  event,  since  the  features  on  the  images  do,  by  inspection, 
have  a  clear,  preferred,  azimuthal  orientation. 

A  key  point  to  remember  is  that  the  unpolarized  images  would  be  similar  to  the  sum  image  (A  +  B)  shown  in  Fig.  2(d). 
Here,  the  details  of  the  river  and  its  tributaries  and  lakes  are  far  less  evident  than  is  the  case  in  the  original  image  A 
(Fig.  2(a)),  the  difference  images  (A-B)  shown  in  Fig.  2(c)  or  the  second  principal  component  image  shown  in  Fig.  2(h). 
The  first  inverse  principal  component  image  shown  in  Fig.  2(i)  is  probably  the  most  useful  image  for  the  delineation 
of  river  features.  Therefore,  the  use  of  polarization  characteristics  as  a  channel  of  information  is  most  helpful  in 
revealing  new  information. 

While  the  above  example  deals  with  hydrological  information,  a  geological  example  is  shown  in  Fig.  3.  The  El  Richat 
structure  in  Mauritania  has  been  the  subject  of  some  discussion.  It  is  uncertain  whether  this  feature  is  the  result  of  the 
impact  of  a  large  meteor,  as  is  the  case  for  Diablo  Canyon.  Differences  in  intensity  and  in  contrast  across  boundaries 
are  clear  from  careful  comparison  of  the  maximum  and  minimum  polarization  images  shown  in  Figs.  3(a)  and  3(b), 
which  constitute  the  polarization  pair.  The  differences  are  clearer  from  reference  to  the  difference  image  (A-B), 
which  is  shown  in  Fig.  3(c),  and  clearer  still  from  reference  to  the  normalized  difference  (i.e.,  polarization)  image 
(A-B)  (A  -r  B),  which  is  shown  in  Fig.  3(d).  These  images  have  been  digitally  enhanced,  and  fine  detail  is  present 
in  the  polarization  image,  which  was  not  present  in  either  components  of  the  original  image  pair.  The  sum  image 
(A  ^  B)  is  shown  in  Fig.  3(e)  and  is  close  to  what  might  be  obtained  by  conventional  (i.e.,  nonpolarized)  photography. 
Comparison  of  the  sum  image  with  the  difference  image,  the  polarization  (normalized  difference)  image  and  the 
ratio  image  shown  in  Fig.  3(f)  clearly  demonstrates  that  the  recording  of  polarization  information  yields  new 
information.  In  particular,  the  relief  topography  is  much  more  evident  in  the  difference,  the  ratio  and  the  polarization 
images  than  is  the  case  in  the  sum  image.  Since  the  illumination  of  shadowed  areas  is  almost  entirely  by  diffuse 
irradiance,  which  is  polarized,  it  is  not  surprising  that  the  upwelling  radiance  from  shadowed  areas  is  strongly  polarized, 
the  degree  of  polarization  depending  upon  the  depth  of  shadow'.  We  have  had  some  positive  indications  that  additional 
information  in  shadowed  areas  may  be  obtained  by  the  digital  analysis  of  polarized  image  pairs  obtained  over  shadowed 
regions.  However,  more  experimental  study  is  needed  before  it  may  be  determined  fully  whether  and  what  extent 
this  is  so. 


5.  CONCLUSIONS 

In  this  study,  we  have  shown  only  a  representative  analysis  performed  on  two  pairs  of  images.  The  polarized  image 
pairs  discussed  were  obtained  over  terrestrial  targets  of  hydrological  and  geological  interest.  We  found  that  the 
differences  between  the  image  pairs  contained  a  significant  amount  of  information  on  contrast,  contrast  boundaries, 
shadow,  and  such  terrestrial  features  as  rivers.  By  examining  the  intensity  image,  which  contained  no  polarization 
information,  we  found  that  polarization  added  greatly  to  the  information  on  surface  features,  such  as  details  on 
an  entire  river  system  or  details  on  topography  that  were  not  clear  in  the  intensity  (unpolarized)  images. 

Combinations  of  the  polarization  images  were  used  as  input  to  a  principal  component  analysis,  which  provided 
some  useful  pictorial  information  not  contained  in  any  of  the  input  images.  Further  benefit  was  obtained  by  performing 
histogram  equalization  on  the  principal  component  images  prior  to  applying  the  inverse  transformation,  thus 
decorrelating  the  input  images.  This  technique  provided  better  definition  of  features  not  observed  in  the  unpolarizcd 
sum  image.  However,  such  approaches  are  empirical  and  cannot,  at  this  time,  be  employed  in  a  deterministic  fashion. 
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By  differencing  the  Fourier  spectra  of  the  polarization  image  pair,  we  found  that  significant  differences  exist  between 
the  images  in  the  population,  azimuthal  distribution  and  spatial  frequency  distribution  of  contrast  boundaries.  If 
contrast  across  a  boundary  is  reduced,  then  there  will  be  a  point  at  which  the  contrast  difference  will  cease  to  be 
observed,  since  it  will  disappear  into  the  noise.  Conversely,  appropriate  use  of  polarization  methods  may  render 
contrast  boundaries  more  readily  observable  by  increasing  contrast  differences  across  the  boundary,  thus  reducing 
clutter  and  improving  signal-to-noise  ratio,  as  observed  by  Coulson,  Whitehead,  and  Campbell.10  In  fact,  for  the 
image  of  the  F.l  Richat  structure  we  found  that  topographic  details  were  considerably  enhanced  by  digitally 
manipulating  the  polarization  image  pairs. 
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Fig.  1.  The  three  states  of  polarization.  Here,  the  mean  electric  vector 
of  recorded  radiance  is  viewed  along  the  direction  of  propagations 
(a)  Shows  a  slate  oj  imear  polarization,  where  the  mean  electric  vector 
does  not  rotate  about  the  direction  of  propagation  with  time,  (b)  Shows 
a  state  of  circular  polarization  in  which  the  ensemble  of  propagated 
electromagnetic  waves  produces  the  same  effect  as  two  orthogona' 
waves  which  are  out  of  phase.  The  amplitude  is  coristant,  but  the  electric 
vector  rotates  with  time  about  the  direction  or  propagation,  (c)  Shows 
the  elliptical  polarization  case  where  the  amplitude  of  the  resultant  wave 
varies  with  time  as  its  plane  rotates  about  the  direction  of  propaga¬ 
tion.  Here,  t,  is  the  orientation  of  the  polarization  ellipse  with  respect 
to  a  reference  direction. 
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I'ig.  2(b).  Minimum  intensity  image  (II)  of  the 
Xiger  River. 
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Fig.  2(g).  First  principal  component  image.  Fig.  2(h).  Second  principal  component  image, 

which  shows  more  detail  of  the  Niger  River 
than  image  A. 


Fig.  2(i).  First  inverse  principal  component 
image  after  histogram  equalization  of  the 
principal  component  images  and  application 
of  the  inverse  transformation  matrix. 
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Fig.  2(j).  Maximum  and  minimum  intensity 
images  constituting  the  polarization  pair  of 
the  Niger  River,  together  with  the  difference 
of  their  power  spectra. 
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Fig.  3(b).  Minimum  intensity  intage  (B)  of  the 
El  Richat  structure. 


Fig.  3(a).  Maximum  intensity  image  (A)  of 
the  El  Richat  structure. 


Fig.  3(c).  Difference  image  (A-B) 


Fig.  3(d).  Degree  of  polarization  image 
(A-B)/ (A  +  B). 


Fig.  3(e).  Sum  image  (A  +  B) 


Fig.  3(f).  Ratio  image  (A/B) 
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